Intriguing magnetism of Fe monolayers on hexagonal transition-metal surfaces 
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Using first-principles calculations, we demonstrate that an Fe monolayer can assume very different 
magnetic phases on hexagonal hep (0001) and fee (111) surfaces of 4d- and 5d-transition metals. Due 
to the substrates' d-band filling, the nearest-neighbor exchange coupling of Fe changes gradually from 
antiferromagnetic (AFM) for Fe films on Tc, Re, Ru and Os to ferromagnetic on Rh, Ir, Pd, and Pt. 
In combination with the topological frustration on the triangular lattice of these surfaces the AFM 
coupling results in a 120° Neel structure for Fe on Re and Ru and an unexpected double-row-wise 
AFM structure on Rh, which is a superposition of a left- and right-rotating 90° spin spiral. 



PACS numbers: 

Triggered by the discovery of the giant- 
magnetoresistance effect and to realize novel spintronic 
device concepts [ij, magnetic nanostructures on surfaces 
have been a focus of experimental and theoretical 
research for more than 20 years now. In particular, 
there has been a tremendous effort to grow ultra-thin 
transition-metal films on metal surfaces and to charac- 
terize and explain their magnetic properties. It is now 
generally believed that these structurally simple systems 
are well understood and more complex nanostructures 
such as atomic chains, clusters, or molecules on surfaces 
have moved into the spotlight of today's research. 

Therefore, it came as a big surprise when it was ex- 
perimentally shown that the prototypical ferromagnet 
Fe becomes a two-dimensional (2D) antiferromagnet on 
the W(001) surface 0. That complex magnetic order 
can be obtained even in single monolayer (ML) mag- 
netic films on non-magnetic substrates has been dramat- 
ically demonstrated by the recent discovery of a spin- 
spiral state for a Mn ML on W(110) [3] and a Mn ML on 
W(001) and a nanoscale magnetic structure for an Fe 
ML on Ir(lll) [5(. Surfaces of Ad- and 5<i-transition met- 
als (TMs) such as W, Re, Ru, or Ir have been particularly 
attractive from an experimental point of view as ultra- 
thin 3d-TM films can often be grown pseudomorphically 
and without intermixing 0, H, 0, E3] • However, there 
has been a controversy in the past about reports concern- 
ing dead magnetic layers and absence of magnetic order 
in ultra-thin films on these surfaces 0,0. The fundamen- 
tal key to many unresolved puzzles may be the itinerant 
character of TMs resulting in competing exchange inter- 
actions beyond nearest-neighbors and higher-order spin 
interactions beyond the Heisenberg model. The latter in- 
teractions have been proposed to play a role in transition- 
metals, however, to our knowledge no unambiguous proof 
of their importance has been given. 

Here, we use first-principles calculations to demon- 
strate that an Fe ML can assume very different magnetic 
phases on a triangular lattice provided by hep (0001) and 
fee (111) surfaces of Ad- and 5<i-transition metals, which 



are also experimentally accessible, ej*. Fe/Ir(lll) 0,0, 
Fe/Ru(0001) 0, or Fe/Pt(lll) [IcJ]. We show that 
the nearest-neighbor exchange interaction, J\, in the Fe 
ML changes continuously from antiferro- to ferromag- 
netic (FM) with filling of the substrate d-band. Due to 
topological frustration on a triangular lattice, AFM cou- 
pling for Fe on Re(0001) and Ru(0001) leads to a Neel 
ground state with angles of 120° between adjacent spins. 
This finding can explain unresolved experimental stud- 
ies reporting the absence of a ferromagnetic signal for 
Fe/Ru(0001) @, 0, [HL[l2|. 

For an Fe ML on substrates such as Ru, Os, Rh, or Ir, 
Ji is small and interactions beyond nearest-neighbors or 
higher-order spin interactions can be relevant. We exem- 
plify this by studying so-called multi-Q states, a super- 
position of symmetry equivalent spin spirals, which are 
degenerate in the Heisenberg model but can gain energy, 
e.g., due to the presence of biquadratic or four-spin inter- 
actions In the Ad-TM substrates, significant mag- 
netic moments are induced by the Fe ML due to their 
high susceptibility. For Fe/Rh(lll), this is a crucial ef- 
fect which stabilizes an unexpected collinear ground-state 
of ferromagnetic double-rows coupling antiferromagneti- 
cally along the [112]-direction, a 2D analog of the antifer- 
romagnetic bilayer state in Fe films on Cu(001) [14{. For 
5d-TM substrates such as Ir or Re, the large spin-orbit 
coupling may cause a significant Dzyaloshinskii-Moriya 
interaction in the Fe ML which can crucially affect the 
magnetic order 0, 0, making these systems a unique 
playground for complex magnetism. 

We have determined the electronic and magnetic prop- 
erties of 1 ML Fe on the hexagonal hep (0001) and 
fee (111) surfaces of Ad- and isoelectronic 5<i-transition 
metals based on density-functional theory. Calculations 
have been carried out in the generalized gradient ap- 
proximation (GGA) to the exchange-correlation func- 
tional [15j using the full-potential linearized augmented 
plane wave (FLAPW) method, as implemented in the 
fleur code [13] • 

The collinear magnetic states were investigated in sys- 
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terns modelled by 6 or 7 layers of 4c?- or 5d-TM substrate 
with hep or fee stacking covered by a pseudomorphic Fe 
monolayer on each side of the films. We have used the ex- 
perimental lattice constants which are very close to the 
values obtained by GGA. The structural relaxation of 
the Fe overlayer has been performed for both fee and hep 
stacking. The non-collinear magnetic states have been 
studied employing an asymmetric film consisting of four 
substrate layers and an Fe monolayer on one side of the 
film at the distance optimized for the collinear (FM or 
AFM) state of lowest energy. For Fe on Rh(lll) we have 
used 6 substrate layers; we found that adding two layers 
of substrate did not influence the spin-spiral dispersion by 
more than 6 meV. The spin spirals have been calculated 
exploiting the generalized Bloch theorem [3]. We have 
used about 100 basis functions per atom for all calcula- 
tions and at least 676 ky points in the two-dimensional 
Brillouin zone (2D-BZ) for the spin-spiral calculations, 
48 k|| points in one quarter of the 2D-BZ for the uudd 
configurations and 32 ky points in the 2D-BZ for the 3Q- 
state requiring a surface unit-cell comprising 4 atoms. 

In order to find the magnetic ground state, we start 
by evaluating the total-energy difference between the FM 
and the row-wise AFM (RW-AFM) configuration, Fig.[Q 
considering hep and fee stackings of the monolayer. Only 
for substrates at the end of the TM-series, Pd and Pt, the 
Fe monolayer prefers fee stacking. On all other substrates 
Fe prefers an hep stacking. Only on Rh and Ir the en- 
ergy difference between fee and hep stacking is sufficiently 
small (9.0 and 7.6 meV/Fe-atom, respectively) to suggest 
the experimental observation of both types [a, [l9[ . 

From Fig. [T] it can be concluded that Fe on substrates 
from the center of the TM-series, Tc and Re, exhibits a 
clear antiferromagnetic behavior and the RW-AFM state 
has the lowest energy. Fe on substrates at the end of 
the TM-series is ferromagnetic. In between we observe 
a gradual change from a strongly AFM behavior to a 
strongly FM one as function of the electron filling of the 
substrate. It is argued that this change in the magnetic 
coupling results from the 3d— Ad and 3d— 5d hybridization 
between the Fe ML and the substrate, which is altered 
by the d-band filling. This argument is supported by the 
fact (i) that the role of the hybridization is also appar- 
ent from the monotonous variation of the Fe magnetic 
moments as one moves through the TM-series and (ii) 
the gradual change from FM to AFM coupling cannot 
be explained on the basis of the changing in-plane lattice 
constant as the comparison with unsupported MLs on 
the respective lattice constants shows a rather stable fer- 
romagnetic value of about 160 meV/Fe-atom, c.f. Fig. Q] 

From the above results we can conclude that the anti- 
ferromagnetic exchange interaction is strong and impor- 
tant for Fe on most of these substrates except Pd and 
Pt on which Fe is clearly ferromagnetic. The antiferro- 
magnetic interaction on a triangular lattice leads to the 
frustration of magnetic interactions and is the origin of 
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FIG. 1: (color online) Total-energy difference between the 
FM (AE > 0) and RW-AFM (AE < 0) configuration for 
Fe MLs on hep (0001) and fee (111) surfaces of 4d- and 5d- 
TMs. Closed and open symbols indicate a favorable hep or 
fee stacking of the Fe ML, respectively. The magnetic mo- 
ment of the Fe atoms, mF C , is given. A_Bapm-fm is nearly 
constant for an unsupported hexagonal Fe ML (UML) with 
the corresponding in- plane lattice constants (dashed line). 

complex magnetic states. Fe on Re or Tc exhibits strong 
antiferromagnetic interactions as shown by the large en- 
ergy gain when assuming a RW-AFM state and the true 
ground state could be a 120°-Neel state. Due to the small 
energy difference between the FM and RW-AFM order 
for Fe on Os, Ru, Rh and Ir many magnetic states have 
to be considered as possible ground states. Therefore, 
we first focus on Fe MLs on Ru(0001) and on Rh(lll) 
as model systems of complex magnetism on a triangular 
lattice [H. 

We study non-collinear magnetic structures by per- 
forming spin-spiral calculations. Flat spin spirals are 
the general solution of the classical Heisenberg model on 
a periodic lattice H — — Y2i<j Jij M» ' > wnere the 
exchange constants determine the strength and the 
type of coupling between local moments at sites i and 
j pointing along the unit vectors Mj and IVLj, respec- 
tively. Spin spirals are characterized by a wave vector 
q and the moment of an atom at site Ri is given by 
Mi(Rj) = M ( cos(q • Rj), sin(q • Rj), 0) , where M is the 
spin moment per atom. By considering spin spirals along 
the high symmetry lines of the 2D-BZ we cover an im- 
portant part of the magnetic phase space. At high sym- 
metry points, we find well-known magnetic states such 
as the FM state at the T-point, the RW-AFM state at 
the M-point, and the 120° Neel state at the K point, 
c.f. Fig. (Ha). 

The calculated total-energy dispersion, E(q), of spin 
spirals for Fe/Ru(0001), Fig. 12(b), shows that the Neel 
state, i.e. for q = K, is the energetically most favorable 
magnetic configuration. The energy gain is 58 meV/Fe- 
atom with respect to the FM state and only 13 meV/Fe- 
atom with respect to the RW-AFM state. The magnetic 
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FIG. 2: (color online) (a) Neel state and uudd-states obtained 
from a superposition of two 90° spin spirals running either 
along TK or TM. (b) Total energy of spin spirals for one hep 
monolayer Fe on Ru(0001) (filled circles) along the high sym- 
metry directions of the 2D-BZ (see inset). Solid line denotes 
a fit to the Heisenberg model up to J5. The 3Q- (down- 
ward triangle) and wudd-state (upward triangle) are included 
for comparison. Magnetic moments of (c) Fe atoms and (d) 
Ru interface atoms, Ru(I) (filled squares) and sub-interface 
atoms, Ru(I-l) (open squares). 

moment in the Fe ML, Fig. (He), depends only weakly on 
q. Interestingly, the moments of the surface and subsur- 
face Ru layer are of similar magnitude and change their 
alignement to the Fe moments, Fig. [2^c)-(d). 

To interpret our calculations, we have mapped the re- 
sults onto a 2D Heisenberg model which allows us to 
determine effective exchange constants Jij between mo- 
ments of the Fe atoms in the ML. Fig. HJb) shows that a 
good fit is obtained by including up to five nearest neigh- 
bors. From the Jy's, given in table HI we conclude that 
the main effect of the substrate is to change the nearest- 
neighbor exchange coupling, Ji, from FM (Ji > 0) for 
Rh and Ir to AFM (J x < 0) for Ru and Re. 

For itinerant magnets such as iron, it is not a priori 
clear that the Heisenberg model, which relies on localized 
magnetic moments, can provide a good description. The 
next higher-order terms beyond the Heisenberg model 
are the biquadratic and four-spin interactions. These 
terms can lead to energy contributions on the order of 
15 meV/Fe-atom 13j which is similar to the energy differ- 



TABLE I: Heisenberg exchange constants for the hep Fe ML 
on different substrates obtained by fitting the total-energy 
dispersion along T — K — M and higher order Heisenberg terms. 
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ence between Neel and RW-AFM state of Fe/Ru(0001). 

In order to find the magnitude of these interactions, 
we consider linear combinations of spin spirals, the multi- 
Q states. The degeneracy of single- and multi-Q states 
within the Heisenberg model is lifted by higher-order 
interactions. By calculating the energy differences be- 
tween suitable single-Q and multi-Q states, we can ob- 
tain values for the nearest-neighbor biquadratic interac- 
tion, B, and four-spin interaction, K . One pair of such 
states are spin spirals at the M-points of the 2D-BZ and 
the 3Q-state constructed from the three independent M- 
points [HI . As a second pair we choose the spin spiral at 
Q3K/4 an( i a superposition of two such spirals with oppo- 
site rotation sense, a so-called uudd-st&te (middle panel 
of Fig. [2ja)). The total-energy differences are 
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marked by squares in Fig. [2^b) and the obtained con- 
stants B and K are given in table U 

The large energy differences between the multi-Q 
and single-Q states of 7 to 11 meV/Fe-atom found for 
Fe/Ru(0001) clearly demonstrate its itinerant character. 
The extracted values for the biquadratic and four-spin 
interaction are of the order of the exchange constants be- 
yond nearest neighbors and cannot be neglected in find- 
ing the ground state. Higher-order interactions play a 
similar role for other substrates, as shown in table fl] 

The Fe ML also induces considerable magnetic mo- 
ments in the Ru substrate due to its high susceptibil- 
ity. The induced magnetic moments of Ru depend on 
the magnetic structure in the Fe ML and their size con- 
tributes to the total energy of the system. This effect 
can be dramatic as seen by comparing the uudd-st&te 
along r — M and Y — K, c.f. right panel of Fig. Ela), to 
the corresponding single-Q states. The energy difference 
due to biquadratic and four-spin interactions should be 
the same in both cases, which is clearly not the case. 
This finding can be traced back to the difference in the 
induced magnetic moments in the substrate. 

We now turn to Fe on Rh(lll) for which the FM state 
was slightly favorable among the collinear states consid- 
ered in Fig. [TJ Our spin-spiral calculations, however, in- 
dicate a non-collinear ground state with a spiral vector q 
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FIG. 3: (color online) (a) Total energy of one hep mono- 
layer Fe on Rh(lll) for spin-spirals along high symmetry lines 
(c.f. Fig. [2J) and multi-Q states. Solid line denotes a fit to 
the Heisenberg model up to J5. Magnetic moments of the 
Fe atoms and the Rh interface atoms (Rh(I), open symbols) 
and sub-interface atoms (Rh(Fl), filled symbols) are shown 
in panels (b) and (c), respectively. 

of about 0.225 x ?f along T - K, about 2.5 meV/Fe- 
atom below the FM state (Fig. [3]). From the fitting 
to the Heisenberg model, c.f. table []] we indeed find a 
very small FM nearest-neighbor exchange constant while 
second- and third-nearest neighbors prefer an AFM align- 
ment of a magnitude comparable to the FM exchange. 

Including the multi-Q states in our search for the 
ground-state results in the uitcM-solution (r — M) as the 
most favorable configuration with a large energy gain of 
29.0 meV/Fe-atom with respect to the single-Q state. 
However, the real magnetic ground-state might be even 
more complex due to the competing interactions involved 
making Fe/Rh(lll) a truly intriguing system and a chal- 
lenge for an experimental investigation. 

In the phase diagram of the 2D Heisenberg model, 
shown in Fig. |4j we can provide a complete picture of 
the substrates' impact on the Fe exchange coupling by in- 
cluding spin spiral calculations for an Fe ML on Re(0001), 
Ir(lll), and Ag(lll). In the J\-J% plane of the diagram, 
Fig. HJa), we see that the d-band filling of the substrate 
drives the system along the line J2 ~ from a Neel con- 
figuration on Re and Ru to the FM solution on Ir and Rh. 
For small J\ , we need to consider also the phase diagrams 
in the J2-J3 plane showing the spin spiral minimum of 
Fe on Rh(lll) in the FKM- direct ion, Fig. H[c), and the 
120° Neel state of Re(0001), Fig. @Jb). In contrast, the 
FM state is the most favorable single-Q state for the hep 
Fe ML on Ir(lll) [2(J. 

In conclusion, we have proposed Fe MLs on hexago- 
nal surfaces of late Ad- and 5<i-TMs as promising sys- 
tems to study experimentally the magnetic interactions 



o 

£ 
'E 

13 



(a) 5 



2. 2. 
=3 


i a: <: 

til en fii ■ ^> ■ 

iS^if * 5 FM 

— 3 


RW-AFM / 


ss r-M 





J, (arb. units) 




FIG. 4: (color online) Phase diagrams of the classical Heisen- 
berg model for a 2D hexagonal lattice for states on the high 
symmetry lines, (a) J1-J2 plane for J3 = 0. J2-J3 plane for 
(b) Ji > and (c) Ji < 0. (Filled symbols denote values ob- 
tained from fits for the Fe ML on 4c?- and 5d-TM substrates.) 

in TMs, e.g. by proving the influence of spin interactions 
beyond the Heisenberg model. Alloying of the substrate, 
e.g. PtRu(0001) [22j], may allow an additional fine tun- 
ing of the degree of disorder in 2D systems as proposed 
in Ref. Q. 
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